Protecting entanglement via the quantum Zeno effect 
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We study the exact entanglement dynamics of two atoms in a lossy resonator. Besides discussing 
the steady-state entanglement, we show that in the strong coupling regime the system-reservoir 
correlations induce entanglement revivals and oscillations and propose a strategy to fight against 
the deterioration of the entanglement using the quantum Zeno effect. 
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The description of decoherence for bipartite entangled 
systems has recently reached notable theoretical [l[ and 
experimental Q results due to the introduction of the 
concept of entanglement sudden death. This describes 
the finite-time destruction of quantum correlations due 
to the detrimental action of independent environments 
coupled to the two subsystems. On the other hand, it 
is well known 0, Q that the interaction with a common 
environment leads to the existence of a highly entangled 
long-living decoherence-free (or sub-radiant) state. At 
the same time, another entangled state exists (orthogo- 
nal to previous one and called super-radiant) that looses 
its coherence faster. In this paper, we discuss how to pre- 
serve this second entangled state without affecting the 
first one by exploiting the quantum Zeno effect in order 
to achieve a complete entanglement survival. 

Specifically, we address an exactly solvable non- 
Markovian model describing two two-level atoms (qubits) 
resonantly coupled to a lossy resonator, which we treat 
through the pseudo-mode approach [1]. As we discuss 
below, this describes, e^., atoms or ions trapped in an 
electromagnetic cavity Q or circuit-QED setups 0, H, Q 
and our results are directly verifiable in both systems. 

We obtain the exact entanglement dynamics as a func- 
tion of the environment correlation time and discuss its 
stationary value, which turns out to be maximal for a 
factorized initial state of the two atoms. In the past, 
the environment induced entanglement generation has 
been discussed in the Born-Markov limit [l^] or for a 
pure dephasing case 



of cavity QED and with superconducting Josephson cir- 
cuits. 

We consider a two-qubits system interacting with a 
common zero-temperature bosonic reservoir. The micro- 
scopic Hamiltonian of the system plus reservoir, is given 
hy H = Hq + _ff int ; with 
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Here, bk is the annihilation operator of quanta of the 
environment, while trji and Uj are the inversion opera- 
tor and transition frequency of the j-th qubit, j ~ 1,2, 
whose interaction with the reservoir is measured by the 
dimensionless constant Oj. This depends on the value 
of the cavity field at the qubit position and can be ef- 
fectively manipulated, e.g., by means of dc Stark shifts 
tuning the atomic transition in and out of resonance. For 
the following discussion, it will prove useful to introduce 
a collective coupling constant ar = {a\ + oi^Y/"^ and the 
relative strengths Vj = Uj/aT (as rf+r'^ ~ 1; '^c take 
only ri as independent). By varying ax, we will explore 
both the weak and the strong coupling regimes. 
For an initial state of the form 



|V(0)) - Coi|l)JO)2 + Co2|0)Jl)2 
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11 1 Here we extend these results ^^"^ time evolution of the total system is given by 

\^it)) = Ci(i)|l)i|0)2|0)B + C2(t)|0)i|l)2|0)B 
+ 5^C,(t)|0)i|0)2|U)£, 



to a dissipative coupling with the environment outside 
the Markovian regime, both for weak and strong cou- 
plings, corresponding to the bad and good cavity limits. 
In particular, in the latter regime, the long memory of the 
reservoir induces entanglement revivals and oscillations. 

Furthermore, we describe a measurement induced 



(4) 



quantum Zeno effect 12] for the entanglement, showing 
that the quite simple procedure of monitoring the pop- 
ulation of the cavity mode leads to a protection of the 
entanglement well beyond its natural decay time. This 
effect too can be tested with slight modifications of al- 
ready existing experimental set-ups, both in the realm 



|lfc)£; being the state of the reservoir with only one ex- 
citation in the k-th mode. Setting J^"'' = uij — ujk, the 
equations for the probability amplitudes take the form 



Ck{t) 
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^*5fc [aie-^'^"*ci(t) + a2e-^^^"*C2(i) . (6) 
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Formally integrating Eq. ([6]) and inserting its solution 
into Eqs. ([5]), one obtains two integro-differential equa- 
tions for the amplitudes Ci 2(0- the continuum limit 
for the environment, and with the introduction of the 
correlation function /(r), and of its Fourier transform 
J{uj) (which is the reservoir spectral density), these two 
equations become 

ci = - / dtif{t - ti) [ajciiti] + aia2C2{ti)] , (7) 
Jo 

C2 = - f dtj{t - h) [alc^ih) + a^a^ciih)] . (8) 

Before discussing the general time evolution, we notice 
that a constant solution can be found, independently of 
the form of the spectral density. Namely, a sub-radiant, 
decoherence-free state exists, that does not decay in time 



r2\l) 
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In the following, wc consider the case in which the two 
atoms have the same Bohr frequency, i.e., uji = uj2 ~ i^o, 
and interact resonantly with a reservoir with Lorentzian 
spectral density 



J{lo) 
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This is, e.g., the case of two atoms interacting with a 
cavity field in presence of cavity losses. Due to the non 
perfect reflectivity of the cavity mirrors, the spectrum 
of the cavity field displays a Lorentzian broadening. In 
this case the correlation function decays exponentially 
/(r) = W'^e^^'^, the quantity 1/A being the reservoir 
correlation time. The ideal cavity limit is obtained for 
A — > 0; in this case one has 



lim J(lu) = W'^Siuj -Luo), 



(11) 



corresponding to a constant correlation function /(t) = 
W^. The system, then, reduces to a two-atom Jaynes- 
Cummings model [13| with a vacuum Rabi frequency 
TZ = arW. On the other hand, for small correlation 
times (with A much larger than any other frequency 
scale), we obtain the Markovian regime characterized by 
a decay rate 7 = 2R? /\. For generic parameter values, 
the model interpolates between these two limits. 

As IV'-) docs not evolve in time, the only relevant time 
evolution is that of its orthogonal, super-radiant, state 



1^^+) -ri|l)i|0)2+r2|0)i|l)2. 
One obtains |V'+(i)> = ^(OI'/'+(0)>: with 
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FIG. 1: (Color online) Stationary concurrence as a function of 
the relative coupling constant ri and of the initial separability 
s of the state, for (a) = 0, and (b) tj) = tt. In the first case, 
the maximum of Cs is achieved for asymmetrical couplings: 
at ri ~ 0.87 for s = 1, and at ri ~ 0.5 for s = —1. In the 
second case the maximum is achieved in correspondence of 



In the {|1), |0)} basis, the reduced density operator for 
the two qubits is given by 
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(10) where, letting /3± = {t/j± \ ^(0)), one has 



ci{t) = r2(3-+n£{t)(3+, 
C2(t) = -nP_+r2£it)P+. 



(14) 
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(16) 



The solution is exact as we have not performed neither 
the Born nor the Markov approximation. We now use 
this result to obtain the dynamics of the qubit entangle- 
ment as measured by the concurrence C (t) [l3| , ranging 
from (for separable states) to I (for maximally entan- 
gled ones). For the density matrix given by Eq. ((M]), the 
concurrence is 



c{t) = 2\cim{t)\ 



(17) 



We begin by noticing that there 
stationary value of C due to the 
the decoherence-free state: Cs = 
Cil^b-)) |(V'-|^^(0))|'=2|rir2| . 
the time evolution of the concurrence as a function of the 
initial amount of entanglement stored in the system, we 
consider initial states of the form (13]) with 



exist a non-zero 
entanglement of 
C(< ^ 00) = 
To better discuss 
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,with- 1 < s < 1. 



where il = y/y^4^. 



Here, the separability parameter s is related to the ini- 
tial concurrence as = 1 — C(0)^. Fig. [T]-(a) displays 
the stationary concurrence versus ri and s. Due to the 
interaction with the cavity field, initial separable states 
(s — ±1) become entangled. In fact, for = 0, the max- 
imum stationary entanglement C™^^ ~ 0.65 is obtained 
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FIG. 2: (Color online) Time evolution of the concurrence in 
the bad cavity limit {R = 0.1, top plots) and good cavity 
limit {R = 10, bottom plots), with (a) s = 1, and (b) s = 0, 
both with </) = 0, for the cases of i) maximal stationary value, 
ri — 0.87 (solid line), ii) symmetrical coupling ri — l/\/2 
(dot-dashed line), and iii) only one coupled atom ri = 0, 1 
(dotted line). The insets show the initial quadratic behavior 
of the concurrence for R = 0.1 

for initially factorized states, i.e. s = ±1. While the 
details depend on the phase (jj, the qualitative picture is 
generic and essentially independent of 0, apart from the 
isolated case of an initial state coinciding with In 
such a situation all of the entanglement initially encoded 
in the qubits remains there for long times. For positive 
Tj, this occurs for (j) ~ n, see Fig. [T]-(b). 

We now look at the entanglement dynamics in the good 
and bad cavity limits, i.e. for i? ^ 1 and i? ^ 1, respec- 
tively, with R = TZ/X. In Fig. [2] we show the concur- 
rence as a function of r = At in the bad (upper row) 
and good (lower row) cavity limits. Wc compare the dy- 
namics of an initially factorized state (s = 1) with the 
one of an initially maximal entangled state (s = 0) for 
four different values of the coupling parameter ri , namely 
n = 0,1/a/2, 0.87,1. The plots for n = and n = 1 
overlap as they both describe a case in which one of the 
two atoms is effectively decoupled. The value ri = 0.87 
corresponds to the case of optimal stationary entangle- 
ment for the initial state s — 1 with (j) = 0. Finally 
ri = 1/ -y/2 corresponds to the case of equal coupling of 
the two atoms with the reservoir. Other values of ri show 
qualitatively similar behavior. 

For weak couplings and/or bad cavity, R = 0.1, and 
for an initially separable state (s = 1), the concur- 
rence increases monotonically up to its stationary value; 
whereas, for initially entangled states, the concurrence 
first goes to zero before increasing towards Cg- The 
strong coupling/good cavity case i? = 10 is more rich 
and presents entanglement oscillations and revival phe- 



nomena for all the initial atomic states. One can prove 
analytically that for maximally entangled initial states 
(s ~ 0) the revivals have maximum amplitude when only 
one of the two atoms is effectively coupled to the cav- 
ity field, i.e. for ri = 0,1. In this case, indeed, the 
system performs damped oscillations between the states 
14'+} and \tp-} which are equally populated at the be- 
ginning. On the other hand, for an initially factorized 
state, the interaction with the cavity field in the strong 
coupling regime generates a high degree of entanglement. 
Indeed, for i? = 10, at r = f ~ 0.31, C attains the value 
C(f) ~ 0.96, at n ~ 0.92 (for s = 1) or ri ~ 0.4 (for 

These entanglement revivals are a truly new result due 
to the memory depth of the reservoir. Very small revivals 
can occur in the Markovian regime [l5| , and in the non- 
Mar kovian regime with independent reservoirs [lit . 
our case, however, the amount of revived entanglement is 
huge, being comparable to the previous maximum. This 
feature only appears in the strong coupling regime and 
with a non-zero environmental correlation time. The sur- 
prising aspect, here, is that an irreversible process such 
as the spontaneous emission is so strongly counteracted 
by the memory effect of the environment, which not only 
creates entanglement, but also lets it oscillate quite a few 
times before a stationary value is reached. 

If we express the initial state of the qubits as a super- 
position of |V'±)j that is IV'(O)) = /3-IV'-) + 13+\'^+), we 
see that, while part of the initial state will be trapped in 
the sub-radiant state another part will decay fol- 

lowing Eq. (fT3|) . Thus, as discussed above, the amount of 
entanglement that survives depends on the specific state 
(and on the value of the rj). In the following we present 
a method which exploits the quantum Zeno effect to pre- 
serve the initial entanglement independently of the state 
in which it is stored. 

We consider the action of a series of nonselective mea- 
surements on the collective atomic system, performed at 
time intervals T, which have the two following properties: 
i) one of the possible measurement outcomes is the pro- 
jection onto the collective ground state \tpg) = |0)i|0)2, 
and ii) the measurement cannot distinguish between the 
states |1)^ 10)2 and \0)^ \1)2- Any procedure fulfilling 
these two conditions will do the job of preserving the 
entanglement. In particular, one could measure the col- 
lective atomic energy [condition ii) holds in this case since 
the transition frequencies are equal] or, more simply, one 
could monitor the state of the cavity: if a photon is found, 
than the qubits have necessarily decayed into \4'g), while 
if no photon is found, than the excitation still resides on 
the atoms. This can be done both in cavity QED set-ups 
(by sending a probe atom through the cavity that can ab- 
sorb the photon) and with superconducting circuits (by 
sending a short measuring voltage pulse to the resonator, 
similarly to Ref . P ) . 

The measurements described above disentangle the 
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FIG. 3: Time evolution of the concurrence, for s = and 
r — l/\/2, in absence of measurements (solid line) and in 
presence of measurements performed at intervals (a) AT = 
0.1,1,5 (dashed, dotted and dot-dashed lines, respectively) 
with R = 0.1 (weak coupling) and (b) XT = 0.001, 0.005, 0.01 
(dashed, dotted, and dot-dashed lines, respectively) with R — 
10 (strong coupling). 



qubits from the reservoir at each time T. Choos- 
ing T such that (*g|p(r)|*g) < 1, it is straightfor- 
ward to prove that the state IV'-) is miaffected and 
that, at the same time, the decay of is slowed 

down. Its survival probabiUty P\^\t) ~ {^jf-\p(t)\^^) 



in presence of N measurements is given by P^"^ (t) ~ 
|/?+(0)|2 exp [-7^(r)t], where t = NT and with an effec- 
tive decay rate 



log [SjT)' 
T 



(18) 



Notice that, in the limit T ^ and N ^ oo, with a 
finite t = NT , jz{T) and the decay is completely 
suppressed. 

Besides affecting the probability P+{t), the projective 
measurements also modify the time evolution of the en- 
tanglement, whose effective dynamics now depends on 
T. Explicitly, the concurrence at time t = NT, after 
performing N measurements, is given by 



C(^)(t) 



(19) 



In Fig.[3]we compare the dynamics of C(t) in absence and 
in presence of measurements performed at various inter- 
vals T for an initially maximal entangled state. Both 
in the weak and in the strong coupling regimes (left 
and right plots, respectively) the presence of measure- 
ments quenches the decay of the concurrence. Thus, we 
have achieved a quantum Zeno protection of entangle- 
ment from the effect of decoherence. Again, decreasing 
the interval between the measurements, C^-^"^ (t) remains 
closer and closer to its initial value. 

To sum up, we discussed an exactly solvable model 
describing two qubits interacting with a non ideal res- 



onator. We analyzed in detail the stationary entangle- 
ment and obtained the entanglement dynamics both in 
the weak and strong coupling limits, showing that en- 
tanglement revivals can appear due to the finite memory 
of such a complex environment. We also investigated 
the quantum Zeno effect for this system, showing that 
the entanglement can be preserved independently of the 
state in which it is encoded, with the help of repeated 
projective measurements. As anticipated above, our re- 
sults apply both to cavity QED experiments with trapped 
atoms and to the case of superconducting circuits, with 
on-chip qubits and resonator. In the first case, it has al- 
ready been demonstrated that both atoms and ions can 
be confined inside high finesse optical cavities and their 
quantum states can be fully controlled Q . In the second 
case, quantum communication between two Josephson 
qubits has been achieved using a transmission line as a 
cavity 
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